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Abstract 

We give a physical interpretation to the multi-polar Erez-Rozen-Quevedo 
solution of the Einstein Equations in terms of bars. We find that each multi- 
pole correspond to the Newtonian potential of a bar with linear density pro- 
portional to a Legendre Polynomial. We use this fact to find an integral 
representation of the 7 function. These integral representations are used in 
the context of the inverse scattering method to find solutions associated to 
one or more rotating bodies each one with their own multi-polar structure. 
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1 INTRODUCTION 



The adequate description of the gravitation field of an astrophysical object 
has been an important subject in both relativistic and Newtonian gravity 
since their origin. The particular case of the gravity associated to an axially 
symmetric body has played a central role in this theme. In Newtonian the- 
ory the gravitational potential of an axially symmetric body can be always 
represented by its usual expansion in terms of Legendre polynomials (zonal 
harmonics). In general relativity we have that the Einstein equations for an 
axially symmetric spacetime M, 

reduces to 

1p,rr + 1p,r/ T + 1p, zz = 0, (2) 

7 [V] = / r[(^ 2 r - f z )dr + 2^, z dz], (3) 

where the functions ip an d 7 depend only on r and z. 

The simplicity of the Weyl equation is rather deceiving since the po- 
tential ip that obeys the usual Newtonian Laplace equation has a different 
meaning in general relativity. The mono-polar solution of the Einstein equa- 
tions (Schwarzschild's solution) is represented in these coordinates by a bar 
of constant linear density of length 2m, i.e., a 'potential' whose multi-polar 
expansions contains all the multi-polar moments beyond the dipole. Newto- 
nian images are useful in general relativity but need to be used carefully [Q]. 
Erez and Rosen || found that in spheroidal coordinates 

x=(R++R-)/(2<r), y=(R +-R_)/(2a) , (4) 
R+ = ^fr 2 + (z + a) 2 , R_ = yV 2 + (z - a) 2 , (5) 

with a = m and x > 1, — 1 < y < 1, the potential ip outside an isolate 
axially symmetric body can be written as 

oo 

^ = - Yl QkQk(x)P k (y), (6) 

k=0 
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where and Qk are the Legendre polynomials and the Legendre functions 
of the second kind, respectively, and g& are constants. The mono-polar terms 
correspond exactly to the potential representative of the Schwarzschild bar 
(<7o = 1, and q^ = 0, for k 7^ 0). The particular case q — 1, q\ = 0, q 2 = q 
and qk = for k > 2 is the solution usually associated to Erez-Rosen and 
studied in text books ||. 

The function 7 for the general solution (§) is rather formidable and 
was first computed by Quevedo 0. Also the physical meaning of the multi- 
polar moments was clarified in ||. This is a rather important point that can 
be settled after the work of Giirsel || that shows the relation between the 
physically oriented definition of moments due to Thorne J7| and collaborators 
and the more geometrically motivated definition of Geroch || and Hansen 
U that are of simpler computation Also the moments for a charged 

stationary generalization of the Erez-Rosen-Quevedo solution were studied 
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The aim of this paper is to give a physical interpretation to each term 
in the multi-polar expansion (|j) in terms of Newtonian potentials associated 
to bars like the one of the above mentioned Schwarzschild bar. This inter- 
pretation is achieved in terms of potentials associated with bars of length 2m 
with densities proportional to Legendre polynomials. This integral represen- 
tation of the function ip can be used to find a simple integral representation 
of the non-trivial function 7. Furthermore, an integral representation of the 
function F used to generate stationary solutions in the context of the inverse 
scattering method (ISM) [12| [13| can be easily found. An example of solution 
generated by this method is a multi-polar deformed Kerr-NUT solution. We 
also discuss a new solution of the Einstein equations that represents two or 
more Kerr-NUT solutions, each one with their own multi-polar deformations, 
located along an axis of symmetry. 

In Sec. II, combining two identities involving Legendre functions we find 
the physical interpretation of @ in terms of the potential associated to bars. 
We also study an integral representation of the 7 function. In Section. Ill, 
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we give a summary of the main formulas of the ISM and we discuss solitonic 
solutions whose seeds are a field of multi-poles (q = 0, q^ ^ 0,k > 0). 
Particular cases are a "deformed" Kerr solutions and multiple "deformed" 
Kerr solutions . We conclude in Sec. IV with a few remarks. 



2 BARS OF VARIABLE DENSITY AND THE 
7 FUNCTION 

In this section we show the relation between the potential ip defined in (^) 
and the Newtonian potential associated to bars. The first step is the Jeffery 
identity |TJj] that relates Legendre functions of second kind (Q n ) with Ferrers 
functions of the first (P^) and second kind (Q™), also known as associated 
Legendre functions, 

r .xcost + ysint + iz (n - m)\ . . . . 

/ QJ ) cosmtdt = 2ti-. f-Q^(zr)PZ n (cos^) cosfrmp) 

J tv c {n + my. 

(7) 

with 

(x, y, z) = c[(f 2 + 1)2 sin $ cos ip, (f 2 + 1) 2 sin ■& sin if, f cos 8]. (8) 
By doing m = 0, if — u, cos$ = v in (H) and using Neumann formula, 



we find 



I dt 

27rP n (v)Q n (u) = - / P n (t')dt' / — — — — , (10) 

I J -1 J— K A + B cos t + C sm t 

with 



A = uv - 1', 

B = (1 - M 2 )5(l - t> 2 )5 COSV9, 

C = (l-u 2 )^(l-v 2 )hm(p. (11) 
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Hence 

Q n (u)P n (v) = - / nl J 12 

2-/-1 ^(« 2 -l)(l-v 2 ) + (m;-t') 2 

Note that this identity is valid for any complex variables u and v. Now 
identifying (u,v) with the spheroidal coordinates (x,y)„ and using the fact 
that the inverse of M) is: 



r = o^(x 2 - 1)(1 - y 2 ), z = axy, (13) 

we find 

Q.MP M = I f (14) 



2 y_<T y/ r 2 _|_ ^ _ ^2 

Since the Newtonian potential of a bar of length 2a with linear density A 
located symmetrically along the z-axis is 

lr 2 + (z- z') 2 



The density associated to the potential —Q n (x)P n (y) is 



X n (z) = \p n {z/a), (16) 



and the Newtonian mass, 



1 f a 

mn= 2J P n( z / (T ) dz - ( 17 ) 

Note that m n = for n ^ and mo = a = m. In other words the mono-polar 
term is the only one that carries mass, fact that is not surprising. 

A direct computation using the identities of shows that the function 
7 can be written as 



7 = 7o + EE» f dz> f 7 dz"P k {z'/a)P l {z"/a) 



k=0 1=0 



r 2 + (z- z')(z - z") 
A(z' - z") 2 R z ,R z n ' 
(18) 
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where 70 is a constant and R z > = yr 2 + (z — z') 2 , etc. This is a singular inte- 
gral representation, the integration need to be taken as the Cauchy principal 
part, i.e., lim^o.e-^o I-a+ri I-a+ei this allows us to extract the finite part of 
the integral, we are left with a constant of the type ln(0) that can be ab- 
sorbed in the constant 70. In summary, the series (0) can be interpreted 
as the infinite sum of the potentials associated with bars of equal length, 
2a = 2m, and linear densities X n (z) = P n (z/a)/2. 



3 STATIONARY SOLUTIONS 

The vacuum Einstein equations for the stationary axial symmetric metric 
ds 2 = g ab (r, z)dx a dx b + e u (dr 2 + dz 2 ) } (19) 



with (x a ) = (x , x ) = (t,(p), and det(g a b) = —r are equivalent to [O 



v 



(rg, r g- 1 ),r + (rg te g-% = 0, (20) 
Inr - - J r[tr(g, r g- 1 - g, z g~ l )dr + 2tr(g >r g~ 1 )]dz, (21) 



where g = (g ab ) and g r l = (g 1 ) jr , etc. 

Equation (|20|) is an integrable system of equations that is closely related 
to the principal sigma model Techniques to actually find solutions 

of these equations are Backlund transformations and the inverse scattering 
method, also a third method constructed with elements of the previous two is 



the "vesture method" , all these methods are closely related JRy . The general 
metric that represents the nonlinear superposition of several Kerr solutions 
with a Weyl solution can be found by using the inverse scattering method 
PH . We find the N soliton solut ion 

9ab = r^^-ErW^/M) (22) 

s=l k,l=l 

g° ab = dzag(-e 2 ^r 2 e~ 2 ^ (23) 
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N, 



kl 



(k) 



cfe 2 ^\ N[ k) =c[ k \r 2 /fi k )e- 2 ^- F ^ 

N 

n 



N 



e^{c N ) 2 T- N2 ' 2 



IT/** 



,(iV+l) 



(24) 
(25) 

2 detT ab (26) 



k=l 



k, I = 1 



^.-^+[(^ fc -^) 2 + r 2 ] 1/2 . 



(27) 



One of most interesting case is the superposition of one black hole with a field 
of multi-poles, i.e., a two soliton solution (N=2). In this case the constants 
Cq and are related to the usual constants by |L^] 



1/^2 



And 



2 2 , j 2 

cr = m + b 



Also, the constants Wj. are 

W 1 = z 1 + a, 



Wo 



a. 



(28) 
(29) 

(30) 



The symbol that appears in e u is an arbitrary constant, m, a, 6 and 
represent the mass, the angular momentum per unit of mass, the NUT 
parameter, and the position along the z-axis of the black hole, respectively. 
The functions F k are the solutions of the system of equations [ III , 



{rd r - Xd z + 2Xd\)F 
(rd z + Xd r )F = rip Z} 



|A=0 



(31) 
(32) 



see also [13] for an integral 



evaluates along the poles fXk, i-e., F k = F\ x =^ k , 
representation of F. 

In summary, the N-soliton solution (|22"|)-(|2"6"|)(even N) represents the su- 
perposition of a Weyl solution characterized by (V'jT), with N/2 rotating 
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bars (Kerr-NUT solutions) located on z = Zk with masses m^, angular mo- 
menta per unit of mass and NUT charges b^. These constants are related 
by equations like fl28|) and ( p9|) . The case N odd was studied in fl"3| . 

The only integration that is left is the system of equations fl3~l|) with the 



initial condition (32). It is easy to check that the function 

2jfe J-° {z' + R z , + X)R zl 1 ; 

satisfies (pi]) and (|32{) . We have excluded the mono-polar term due to the 
fact that the massive bar will appear as a result of the application of the 
ISM. The last expression can be obtained from the fact that the solution to 
(|3T|)-(p2D associated with the particular solution of Laplace equation 1/R Z > 
is [z + R Z ')/{z + R z > + A). 

Now let us comeback to the two soliton solution. We have that the 
ISM produces a Kerr-NUT bar of length 2a located along the z-axis with its 
center on z — Z\. Really the this bar is independent of the collection of bars 
that describes the multi-polar structure. To have a solution that represent 
an isolated body with a multi-polar structure we need to have bars of the 
same size located on the same position. This is achieved taking z\ = i.e., 
centering the Kerr bar on the origin of the z-axis and instead of taking a = m 
as in the static solution we need a 2 = m 2 + b 2 — a 2 . This solution is a special 



case of the one discussed in |]TTJ that has electric charge. 

The fact that the integral representation is written in cylindrical coor- 
dinates can be used to construct solutions representing two or more bodies 
each one with its independent multi-polar structure located along the z-axis 
in a simple way. In other words these coordinates allow us the separation of 
the different contributions without a considerable amount of work. In prin- 
ciple these configurations are not gravitationally stable, but the nonlinear 
interaction produces struts and/or strings that keep the bodies apart, see for 
instance fl6| . 

To be more specific the metric functions associated to two rotating 
bodies with multi-polar structure can be obtained from applying the ISM to 
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the Weyl solution 

- * - Eb£° S" PJ lW ld >' + ® f ^rVl. < 34 > 

fc=i -Ry i?y 

where 

= J? + (z-z k - z'f- (35) 

Thus ([34]) represents two collections of massless bars each one formed by bars 
of equal size 2a x and 2o 2 located along the z-axis around the points z — z x 
and z = z 2 . The respective function F is, 



(1) r P k (z'/a)(z' + R^) JJ , _ (2) r P k {z'/a){z' + Rf) 



-2F = r 1 r > . 

fc =i J-n (z> + R$ + A)i#> ^ (*' + + A)^ } 

' (36) 

To obtain the two Kerr bars we construct a four soliton solution N = 4, 
in this case we will have eight constants C% that are related to the usual 
constants by [[Lj]] 



C\Cl-ClC\ = o x , C\Cl + ClC{ = m x (37) 

C\C\ — CqCq = —bi, G X C\ + CqCq = a± (38) 

CfC^ - C 3 Ct = a 2 , ClC* + ClCf = m 2 (39) 

ClCt - ClC A Q = -b 2 , C\C\ + C 3 C 4 = a 2 . (40) 



And 



af — ml + bj- a\, o\ = m\ + b\ - a\ (41) 
Also, the constants W k are 

Wx = z x + a u W 2 = z x - at (42) 
W 3 = z 2 + 02, iy 4 = z 2 - a 2 (43) 

The generalization for the case of n = N/2 isolated bodies each one with 
their own multi-polar structure is obvious. 
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4 DISCUSSION 



Multiple Kerr solutions with or without multi-polar structure has been dis- 
cussed by several authors in different contexts and using a variety of solution 
generating techniques fl8 ]. The four soliton representation of two Kerr so- 



lutions each one with their own multipolar deformations presented here has 
a clear meaning. Note that the multi-polar structure of the solution is gen- 
erated by the multi-poles of the Erez- Rosen- Quevedo solution that have a 
precise meaning in general relativity. 

For the single deformed Kerr solution [solution (22)- (27) with N = 2 
and papameters (28)-(30)] as well as for the static case we did not study 
the elimination of conical singularities. The ISM like other similar methods 
introduce singular transformations like t — > t + ap via rotation of tetrads. 
This point together with the study of the horizons, the stability of the two 
body configuration and its related struts and rotating strings will be studied 
in a wider context. The study of struts and strings associated to two Kerr 
metrics can be found in [[17]]. Even though the strut are not physical objects 
(they have negative energy density) they can be used to extract useful infor- 
mation like the rate of emission of gravitational radiation of head-on collision 



of rotating black holes [IS 



When we choose the constants in such a way tha the NUT parame- 
ters bk are zero the ISM, for even N, generates asymptotically flat solutions as 
long as the seed solutions are also asymptotically flat fll3f . This is the case for 
the solutions generated from (|6]). Again in this analysis conical singularities 
are excluded. 
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